A geometric multigrid library for quadtree/octree AMR grids coupled to MPI-AMRVAC

We present an efficient MPI-parallel geometric multigrid library for quadtree (2D) or octree (3D) grids with adaptive refinement. Cartesian 2D/3D and cylindrical 2D geometries are supported, with second-order discretizations for the elliptic operators. Periodic, Dirichlet, and Neumann boundary conditions can be handled, as well as free-space boundary conditions for 3D Poisson problems, for which we use an FFT-based solver on the coarse grid. Scaling results up to 1792 cores are presented. The library can be used to extend adaptive mesh refinement frameworks with an elliptic solver, which we demonstrate by coupling it to MPI-AMRVAC. Several test cases are presented in which the multigrid routines are used to control the divergence of the magnetic field in magnetohydrodynamic simulations

A geometric multigrid library for quadtree/octree AMR grids coupled to MPI-AMRVAC

We present an efficient MPI-parallel geometric multigrid library for quadtree (2D) or octree (3D) grids with adaptive refinement. Cartesian 2D/3D and cylindrical 2D geometries are supported, with second-order discretizations for the elliptic operators. Periodic, Dirichlet, and Neumann boundary conditions can be handled, as well as free-space boundary conditions for 3D Poisson problems, for which we use an FFT-based solver on the coarse grid. Scaling results up to 1792 cores are presented. The library can be used to extend adaptive mesh refinement frameworks with an elliptic solver, which we demonstrate by coupling it to MPI-AMRVAC. Several test cases are presented in which the multigrid routines are used to control the divergence of the magnetic field in magnetohydrodynamic simulations

Radiation-hydrodynamics with MPI-AMRVAC: Flux-limited diffusion

Radiation controls the dynamics and energetics of many astrophysical environments. To capture the coupling between the radiation and matter, however, is often a physically complex and computationally expensive endeavor. Aims. We sought to develop a numerical tool to perform radiation-hydrodynamics simulations in various configurations at an affordable cost. Methods. We built upon the finite volume code MPI-AMRVAC to solve the equations of hydrodynamics on multi-dimensional adaptive meshes and introduce a new module to handle the coupling with radiation. A non-equilibrium, flux-limiting diffusion approximation was used to close the radiation momentum and energy equations. The time-dependent radiation energy equation was then solved within a flexible framework, fully accounting for radiation forces and work terms and further allowing the user to adopt a variety of descriptions for the radiation-matter interaction terms ("opacities"). Results. We validated the radiation module on a set of standard test cases for which different terms of the radiative energy equation predominate. As a preliminary application to a scientific case, we calculated spherically symmetric models of the radiation-driven and optically thick supersonic outflows from massive Wolf-Rayet stars. This also demonstrates our code's flexibility, as the illustrated simulation combines opacities typically used in static stellar structure models with a parametrized form for the enhanced line-opacity expected in supersonic flows. Conclusions. This new module provides a convenient and versatile tool for performing multi-dimensional and high-resolution radiative-hydrodynamics simulations in optically thick environments with the MPI-AMRVAC code. The code is ready to be used for a variety of astrophysical applications, where our first target is set to be multi-dimensional simulations of stellar outflows from Wolf-Rayet stars

Magnetized Accretion-Ejection Structures: 2.5D MHD simulations of continuous Ideal Jet launching from resistive accretion disks

We present numerical magnetohydrodynamic (MHD) simulations of a magnetized accretion disk launching trans-Alfvenic jets. These simulations, performed in a 2.5 dimensional time-dependent polytropic resistive MHD framework, model a resistive accretion disk threaded by an initial vertical magnetic field. The resistivity is only important inside the disk, and is prescribed as eta = alpha_m V_AH exp(-2Z^2/H^2), where V_A stands for Alfven speed, H is the disk scale height and the coefficient alpha_m is smaller than unity. By performing the simulations over several tens of dynamical disk timescales, we show that the launching of a collimated outflow occurs self-consistently and the ejection of matter is continuous and quasi-stationary. These are the first ever simulations of resistive accretion disks launching non-transient ideal MHD jets. Roughly 15% of accreted mass is persistently ejected. This outflow is safely characterized as a jet since the flow becomes super-fastmagnetosonic, well-collimated and reaches a quasi-stationary state. We present a complete illustration and explanation of the `accretion-ejection' mechanism that leads to jet formation from a magnetized accretion disk. In particular, the magnetic torque inside the disk brakes the matter azimuthally and allows for accretion, while it is responsible for an effective magneto-centrifugal acceleration in the jet. As such, the magnetic field channels the disk angular momentum and powers the jet acceleration and collimation. The jet originates from the inner disk region where equipartition between thermal and magnetic forces is achieved. A hollow, super-fastmagnetosonic shell of dense material is the natural outcome of the inwards advection of a primordial field.Comment: ApJ (in press), 32 pages, Higher quality version available at http://www-laog.obs.ujf-grenoble.fr/~fcass

General-relativistic Resistive Magnetohydrodynamics with Robust Primitive-variable Recovery for Accretion Disk Simulations

Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black hole at the center of the galaxy M87 by the Event Horizon Telescope, and the detection of an orbiting "hot spot" nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the development of novel numerical methods to understand the underlying plasma microphysics. Non-thermal emission related to such hot spots is conjectured to originate from plasmoids that form due to magnetic reconnection in thin current layers in the innermost accretion zone. Resistivity plays a crucial role in current sheet formation, magnetic reconnection, and plasmoid growth in black hole accretion disks and jets. We included resistivity in the three-dimensional general-relativistic magnetohydrodynamics (GRMHD) code BHAC and present the implementation of an implicit–explicit scheme to treat the stiff resistive source terms of the GRMHD equations. The algorithm is tested in combination with adaptive mesh refinement to resolve the resistive scales and a constrained transport method to keep the magnetic field solenoidal. Several novel methods for primitive-variable recovery, a key part in relativistic magnetohydrodynamics codes, are presented and compared for accuracy, robustness, and efficiency. We propose a new inversion strategy that allows for resistive-GRMHD simulations of low gas-to-magnetic pressure ratio and highly magnetized regimes as applicable for black hole accretion disks, jets, and neutron-star magnetospheres. We apply the new scheme to study the effect of resistivity on accreting black holes, accounting for dissipative effects as reconnection

Radiatively inefficient MHD accretion-ejection structures

We present magnetohydrodynamic simulations of a resistive accretion disk continuously launching transmagnetosonic, collimated jets. We time-evolve the full set of magnetohydrodynamic equations, but neglect radiative losses in the energetics (radiatively inefficient). Our calculations demonstrate that a jet is self-consistently produced by the interaction of an accretion disk with an open, initially bent large-scale magnetic field. A constant fraction of heated disk material is launched in the inner equipartition disk regions, leading to the formation of a hot corona and a bright collimated, super-fastmagnetosonic jet. We illustrate the complete dynamics of the ``hot'' near steady-state outflow (where thermal pressure $\simeq$ magnetic pressure) by showing force balance, energy budget and current circuits. The evolution to this near stationary state is analyzed in terms of the temporal variation of energy fluxes controlling the energetics of the accretion disk. We find that unlike advection-dominated accretion flow, the energy released by accretion is mainly sent into the jet rather than transformed into disk enthalpy. These magnetized, radiatively inefficient accretion-ejection structures can account for under-luminous thin disks supporting bright fast collimated jets as seen in many systems displaying jets (for instance M87).Comment: Astrophysical Journal (in press). Figures are missing due to file size restrictions. To have the complete paper just click on http://www-laog.obs.ujf-grenoble.fr/~fcasse/MS56638.pd

MPI-AMRVAC: A parallel, grid-adaptive PDE toolkit

We report on the latest additions to our open-source, block-grid adaptive framework MPI-AMRVAC, which is a general toolkit for especially hyperbolic/parabolic partial differential equations (PDEs). Applications traditionally focused on shock-dominated, magnetized plasma dynamics described by either Newtonian or special relativistic (magneto)hydrodynamics, but its versatile design easily extends to different PDE systems. Here, we demonstrate applications covering any-dimensional scalar to system PDEs, with e.g. Korteweg-de Vries solutions generalizing early findings on soliton behaviour, shallow water applications in round or square pools, hydrodynamic convergence tests as well as challenging computational fluid and plasma dynamics applications. The recent addition of a parallel multigrid solver opens up new avenues where also elliptic constraints or stiff source terms play a central role. This is illustrated here by solving several multi-dimensional reaction-diffusion-type equations. We document the minimal requirements for adding a new physics module governed by any nonlinear PDE system, such that it can directly benefit from the code flexibility in combining various temporal and spatial discretisation schemes. Distributed through GitHub, MPI-AMRVAC can be used to perform 1D, 1.5D, 2D, 2.5D or 3D simulations in Cartesian, cylindrical or spherical coordinate systems, using parallel domain-decomposition, or exploiting fully dynamic block quadtree-octree grids

MPI-AMRVAC: A parallel, grid-adaptive PDE toolkit

We report on the latest additions to our open-source, block-grid adaptive framework MPI-AMRVAC, which is a general toolkit for especially hyperbolic/parabolic partial differential equations (PDEs). Applications traditionally focused on shock-dominated, magnetized plasma dynamics described by either Newtonian or special relativistic (magneto)hydrodynamics, but its versatile design easily extends to different PDE systems. Here, we demonstrate applications covering any-dimensional scalar to system PDEs, with e.g. Korteweg–de Vries solutions generalizing early findings on soliton behavior, shallow water applications in round or square pools, hydrodynamic convergence tests as well as challenging computational fluid and plasma dynamics applications. The recent addition of a parallel multigrid solver opens up new avenues where also elliptic constraints or stiff source terms play a central role. This is illustrated here by solving several multi-dimensional reaction–diffusion-type equations. We document the minimal requirements for adding a new physics module governed by any nonlinear PDE system, such that it can directly benefit from the code flexibility in combining various temporal and spatial discretization schemes. Distributed through GitHub, MPI-AMRVAC can be used to perform 1D, 1.5D, 2D, 2.5D or 3D simulations in Cartesian, cylindrical or spherical coordinate systems, using parallel domain-decomposition, or exploiting fully dynamic block quadtree-octree grids